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PACS 64 . 60 . Ht - Dynamic critical phenomena 

Abstract -Based on the Hamiltonian equation of motion of the </> 4 theory with quenched disorder, 
we investigate the depinning phase transition of the domain-wall motion in two-dimensional mag- 
nets. With the short-time dynamic approach, we numerically determine the transition field, and 
the static and dynamic critical exponents. The results show that the fundamental Hamiltonian 
equation of motion belongs to a universality class very different from those effective equations of 
motion. 
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Introduction. — In recent years the dynamics of elas- 
tic systems in disordered media has been a focus of the- 
oretical and experimental studies [THS]. Particularly, the 
magnetic domain-wall dynamics is an important topic in 
magnetic devices, nanomaterial, thin films, and semicon- 
ductor [ThTu]. At zero temperature, the domain- wall mo- 
tion exhibits a depinning phase transition. Due to the 
energy barriers created by the disorder, the domain wall 
is pinned and the velocity of the domain wall remains zero 
up to a transition field H c [TTlfM] . At finite temperature, 
the phase transition is softened. 

At present, most theoretical approaches to the domain- 
wall dynamics in ferromagnetic materials are typically 
based on the Edwards- Wilkinson equation with quenched 
disorder (QEW). This equation is a phenomenological 
model and detailed microscopic structures and interac- 
tions of real materials are not concerned. To understand 
the domain-wall motion at a microscopic level, one may 
build lattice models based on microscopic structures and 
interactions. Very recently, with Monte Carlo simulations, 
the depinning transition and relaxation-to-creep transition 
have been carefully examined in the two-dimensional (2D) 
random-field Ising model with a driving field (DRFIM) 
[15HT8] . The results indicate that the QEW equation and 
DRFIM model are not in a same universality class, mainly 
due to the presence of overhangs and islands in the DR- 
FIM model. Both the QEW equation and Monte Carlo 
dynamics of the DRFIM model are effective theories, al- 
though they may explain the experimental phenomena to 
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some extent. It is desirable to explore the fundamental 
theory of the domain- wall motion. 



In general, it is believed that statistical mechanics is 
originated from the fundamental equations of motion for 
many body systems, even though there has not been a gen- 
eral proof. In statistical mechanics, ensemble theories and 
stochastic equations of motion are effective descriptions of 
static and dynamic properties of the statistical systems, 
respectively. In contrast, whether the fundamental equa- 
tions of motion, e.g., Newton, Hamiltonian and Heisenberg 
equations etc, could really produce the same results of sta- 
tistical mechanics remains open [T9rf23]. Numerical solu- 
tions of the Hamiltonian equations of the (f> 4 theory and 
XY model, show that both equilibrium and nonequilib- 
rium properties of statistical systems could be described 
by the fundamental equation of motion, at least for the 
order-disorder phase transition [23H26] . The Hamiltonian 
equation of the (j) 4 theory falls into a same universality 
class of the Monte Carlo dynamics of the Ising model. 



The purpose of this paper is to examine the depin- 
ning phase transition in two-dimensional magnets, with 
the Hamiltonian equations of motion of the 2D </> 4 the- 
ory with quenched disorder, and to compare the results 
with the effective dynamics such as the QEW equation 
and Monte Carlo dynamics. The short-time dynamic ap- 
proach is adopted to numerically extract the static and 
dynamic exponents. 



p-1 



R. H. Dong, B. Zheng and N. J. Zhou 



The model. — The Hamiltonian of the 2D 4 theory 
with quenched disorder on a square lattice is 

i fi 

+^g^-h i fa-H<f }i ], (i) 

with ni = <pi. The quenched random field hi is uniformly 
distributed within an interval [—A, A], and H is a homo- 
geneous external driving field. In this paper, we fix the 
parameters m 2 = 6, g = 1.8 and A = 6. The Hamiltonian 
equation of motion is then written as 

fa = ]T)(<W + fc-n - 2 ^) + ™?(j>i - + h, + H. (2) 

Since energy is conserved, solutions of the equation in the 
long-time regime are assumed to generate a microcanon- 
ical ensemble. The temperature could be defined as the 
averaged kinetic energy in the long-time regime, but it 
may not be proper in the short-time regime. Instead, the 
energy density could be the control parameter, since it is 
conserved and can be input from the initial state. For 
hi = and H = 0, the ground state of the Hamiltonian in 
Eq. (TTJ) is given by 7Tj = and <fii = Co = -J 6m 2 jg « 4.47. 

We consider a dynamic process with a semiordercd ini- 
tial state with a perfect domain wall in the y direction. 
As time evolves, the domain wall propagates and rough- 
ens, therefore it looks like an interface. Simulations are 
performed on a square lattice, with a linear size 2L in the 
x direction and L in the y direction. Antiperiodic and 
periodic boundary conditions are adopted in the x and 
y directions, respectively. Following Refs [TTJ[T31[^7] , we 
rotate the square lattice such that the initial domain wall 
orients in the (11) direction. This avoids the pinning effect 
irrelevant for the quenched disorder. 

To prepare a semiordercd configuration, we first set the 
initial kinetic energy to zero, i.e., 7Ti(0) = 0. The config- 
uration then consists of two fully ordered domains with 
the same magnitude |</>i(0)| = Co and opposite signs. The 
energy of the semiordered state differs from that of the 
ground state only by the domain wall. In the limit of 
L — > oo, this energy difference is infinitcsimally small. In 
other words, the temperature is close to zero. 

After the initial state is prepared, we simply discretize 
fa by [(j) l (t + At) + (j) l {t-At)-2(j) l (t)}/(At) 2 , and update 
the equation of motion up to t = 10 4 with At = 0.02. Ad- 
ditional simulations with At = 0.05 show that At = 0.02 
is sufficiently small. At least 10000 samples of disorder re- 
alization are used for average. Most simulations are per- 
formed with L = 512, and simulations with different L 
confirm that the finite-size effects are negligible. 

To study the depinning transition of the domain inter- 
face, we first introduce a height function 

1 L 

h(y,t) = z J2^ v (t). (3) 

x=l 



Here <f> xy (t) denotes the field at site (x,y). We may also 
map 4> xy to a Ising spin, i.e., S xy = Sgn{4> xy ), and then 
define the height function as in the DRFIM model jl5lfT8] . 
Our simulations show that after a microscopic time scale 
t mic , these two definitions yield the same results. 

With the height function, the average velocity of the 
interface can be calculated 

where < • • • > includes the statistical average and average 
over y. The roughness function is defined as 

u 2 {t) = (h(y,t) 2 ) - (h(y,t)) 2 . (5) 

Let us assume that the depinning transition is of sec- 
ond order. Based on the dynamic scaling arguments sup- 
ported by the renormalization-group calculations, one ex- 
pects that the order parameter v obeys the scaling form 

v{t, r, L) = b-^"v(b-% ^/"t, b^L). (6) 

Here /3 and v are the static exponents, z is the dynamic 
exponent, and r = (H — H c )/H c . The parameter b is an 
arbitrary rescaling factor and L is the lattice size. In this 
dynamic system, 

m ~ t i/z (?) 

is assumed to be the nonequilibrium correlation length. 
Setting b ~ the dynamic scaling form is rewritten as 

v(t, t, L) = t~^ vz v{l, t 1/yz r, t-V*L). (8) 

In the short-time regime, i.e., the regime with £ (t) <C L, 
the finite-size effect is negligibly small, 

v(t,T)=t-P/ vz v(t 1 /» z T). (9) 

At the transition point r = 0, a power- law behavior is 
expected, 

v (t) = r^ vz . (10) 

With Eq. ©, the transition field H c may be located by 
searching for the best power-law behavior of v(t,r) [HI 
29 , 32| . The critical exponent P/vz is then estimated from 
Eq. (|10[) . To determine the critical exponent 1/vz, we 
compute the logarithmic derivative of v(t,r) [151129] 

d T Zn«(t,T)| T=0 ~^*. (11) 

As the domain wall propagates, it also roughen as an 
interface. According to the definition of the roughness 
function lo 2 (t) in Eq. ([5]) , it includes the dynamic evolution 
of the bulk. This may lead to a deviation from the scaling 
behavior of ui 2 (t) in the short-time regime. Therefore we 
redefine a pure roughness function Duj 2 (t)hy eliminating 
the contribution from the bulk 

Du 2 (t) = u 2 (t)-u 2 b (t), (12) 
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where 0J 2 (t) is actually the line susceptibility of the bulk 
in the x direction. uj 2 (t) can be computed from simula- 
tions starting from the fully ordered state. In fact, uJ 2 (t) 
quickly stabilizes to a constant. Thus, for a sufficiently 
large lattice, the pure roughness function should scale as 

55EBH3HI32I 

Du\t) ~ m 2( , (13) 

where £ is the global roughness exponent. 

To independently estimate the exponent z, we introduce 



F(t) = [M {2) (t) - M(t) 2 ]/uj 2 (t). 



(14) 



Here M(t) is the global magnetization and M^ 2 \t) is its 
second moment. In fact, F(t) is nothing but the ratio 
of the planar susceptibility and line susceptibility. Since 
ui 2 (t) describes the fluctuation in the x direction and 
M^ 2 \t) —M(t) 2 includes those in both x and y directions, 
the dynamic scaling behavior of F(t) should be [TBI [29] 



F(t) ~ m/L ~ t x '*/L. 



(15) 



A more informative quantity is the height correlation 
function 



C(r,t) = ([h(y + r,t)-h(y,t)} 2 ). 



(16) 



Similar to the roughness function, one needs to eliminate 
the evolution of the bulk, Cb{r,t). The pure height cor- 
relation function DC(r, t) = C(r, t) — C&(r, t) should obey 
the standard scaling form 



DC(r,t) 



£2(C-Cl<»e)A r 2Clo 



r « m 

m < r 



(17) 



with Q oc being the local roughness exponent. 



Numerical solutions. — The Hamiltonian equation 
of motion is deterministic in character, very different from 
the stochastic dynamics such as the QEW equation and 
Monte Carlo simulation. We should first verify the exis- 
tence of the pinning state. In Fig. [TJ the roughness func- 
tion co 2 (t) and velocity v(t) of the domain interface are 
plotted as functions of t for the disorder strength A = 6 
with a weak driving field H = 0.1. Obviously, uj 2 ap- 
proaches a constant and v drops down to zero within a 
certain period of time. These results clearly indicate that 
the domain wall is pinned by the disorder for a weak driv- 
ing field. With a strong driving field H, the velocity will 
reach a constant, and the domain wall is depinncd. Inbe- 
tween these two states, there may exist a depinning phase 
transition. 

In our simulations, three disorder strengths, A = 2, 
4 and 6 have been examined. However, for the weaker 
disorder, the depinning transition is not typically of second 
order, it is likely of first order or in the crossover regime. 
The power-law behavior in Eq. (fTOf can not be identified. 
This is similar to the DRFIM model with a weaker disorder 
[15 16] . Therefore, we concentrate our attention on A = 6, 
at which the depinning transition is of second order. 



In order to determine the transition field H c , we perform 
simulations with different values of H to search for the best 
power-law behavior of the velocity v(t, r). In Fig. [2] v(t, r) 
is plotted as a function of t for different driving field H 
on a log-log scale. A rather good power-law behavior is 
observed around H = 0.983. Careful analysis of the data 
locates the transition field H c = 0.9828(3). Both curves 
of L = 512 and 1024 nicely stabilize to a same power-law 
behavior, suggesting that the finite-size effect is negligibly 
small. The exponent (3/vz = 0.773(3) is measured from 
the slope of the curve at H c , according to Eq. (fT0"j) . 

To extract the dynamic exponent z, F(t) defined in 
Eq. (fT4|) is displayed at H c on a log-log scale in Fig. [3] 
A power-law growth of F(t) is observed from relatively 
early times. Based on Eq. (|15[) . a direct measurement of 
the slope of the curve gives 1/z — 0.513(4). 

To estimate the exponent v, we compute approximately 
the logarithmic derivative of the velocity in Eq. (fTTj) . In 
Fig. |H d T \nv(t,T) is plotted at H c on a log-log scale. A 
power-law behavior is detected at later times, and a direct 
measurement from the slope of the curve yields 1/vz = 
0.711(7). Taking into account the corrections to scaling at 
early times, we fit the curve to the form of d T In v(t, r) ~ 
t Y l vz {\ + c/t). As shown in Fig.[4] the fitting to the curve 
is good, and it refines the measurement of the exponent 
to 1/vz = 0.731(9). 

From the measurements of ft/vz, 1/vz, and 1/z, we 
compute the critical exponents (3 = 1.06(1), v = 0.702(11) 
and z = 1.95(2). 

From Eqs. (j!3[) and (fT5|) . one expects the scaling form 
Du 2 (t) ~ F(t) 2 ^. Our data analysis shows that deter- 
mination of £ from this scaling behavior is slightly bet- 
ter than from Eq. ([13)) . In Fig. O therefore, the pure 
roughness function Duj 2 (t) is plotted as a function of F(t) 
at H c on a log-log scale. A power-law behavior is ob- 
served after a transient regime. From the slope of the 
curve one measures £ = 0.849(3). To further refine the 
result, we consider a power-law correction in the form of 
Du 2 {t) ~ F{t) 2 i(l + c/F(t)), it leads to C = 0.865(5). 
Compared with Dui 2 (t), u> 2 (t) obviously suffers from a 
stronger correction to scaling at early times. 

In Fig. [HI the pure height correlation function DC (r, t) 
is displayed for different r at H c . For a large r ^> £(f), 
e.g., r = 256, one extracts the exponent 2£/z = 0.885(6) 
by Eq. dTTJ), consistent with ( = 0.865(5) from Fig. [5] For 
a small r DC(r,t) should be independent of t if 

C = Cloc, according to Eq. JIT]). In Fig. [BJ DC(r,t) at 
r = 2 clearly increases with t. A power-law behavior is 
observed with an exponent 2(( — Qoc)/z — 0.196(3). One 
thus obtains (i oc = 0.674(6). Additionally, one may also 
estimate (i oc from DC(r,t) ~ r 2 ^ loc at a large t, e.g., t = 
10 4 . The result 2( loc = 1.35(1) agrees with ( loc = 0.674(6) 
from Fig. [6] 

In Table [lj all the critical exponents measured for the 
Hamiltonian equation of motion of the 2D (j> 4 theory with 
quenched disorder are compared with those for the Monte 
Carlo dynamics of the 2D DRFIM model and the 2D QEW 
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equation. Obviously, the <j> theory belongs to a univer- 
sality class very different from that of either the DRFIM 
model or the QEW equation. The exponent /3 of the </> 4 
theory is much larger. By the definition in the steady 
state, the velocity v(t) ~ t 13 . In the regime H > H c , 
therefore, the external field easily drives the domain wall 
propagating. The local roughness exponent Q oc of the </> 4 
theory is smaller than 1, similar to the DRFIM model, 
and different from the QEW equation. In the latter case, 
Cioc ~ 1, since the domain wall is simply assumed to be 
single-valued. However, it is somewhat surprising that 
the global roughness exponent £ of the 4 theory is also 
smaller than 1, different from both the DRFIM model and 
QEW equation. In the DRFIM model, it is believed that 
overhangs and islands lead to £ > 1. For the QEW equa- 
tion, it is a kind of inconsistence that ( > 1 contradicts 
the single- valued domain wall. 

To understand the roughening phenomena of the </> 4 the- 
ory with quenched disorder, the snapshot of the domain 
interface is displayed at H c in Fig. [7] As time evolves, 
the domain wall propagates and roughens. Overhangs are 
created, but not so prominent as in the DRFIM model 
fT5"] . Especially, islands can hardly be seen, except for a 
number of small ones in the early times. Based on these 
facts, one may argue that the global roughness exponent 
( of the 4> A theory should be smaller than that of the DR- 
FIM model. However, it remains somewhat puzzling why 
it is smaller than 1. 

Experiments of the domain interface exactly at zero 
temperature do not exist. For T > and < H < H c , it 
is reported that the local roughness exponent Q oc — 0.7(1) 
and 0.69(7) in the experiments with ultrathin Pt/Co/Pt 
films [2JEJ, and Qoc = 0.78(1) with Co 28 Pt 72 alloy films 
[33] . Our numerical value Q oc = 0.674(6) is compatible 
with these experimental results. We do not find direct 
experimental measurements of the global roughness ex- 
ponent £ of the domain interface, but other experiments 
such as those on kinetic surface roughening report results 
of both C > 1 and C < 1 HO- 

Finally we remind ourself that at the order-disorder 
phase transition, the Hamiltonian equation of motion of 
the 4> 4 theory falls into a same universality class of the 
Monte Carlo dynamics of the Ising model [25] . Very proba- 
bly this is because the order-disorder phase transition is an 
equilibrium phase transition, and only relaxation dynam- 
ics is concerned. However, the depinning phase transition 
of the domain-wall motion is a dynamic phase transition. 
It is usually believed that the Monte Carlo dynamics may 
not quantitatively describe the dynamic transportation or 
propagation. Therefore the Monte Carlo dynamics fails to 
produce the results of the Hamiltonian equation of motion 
at the depinning phase transition. 

Conclusion. — We have numerically solved the 
Hamiltonian equation of motion of the 2D (j) theory 
with quenched disorder, and carefully examined the depin- 
ning phase transition of the domain-wall motion in two- 



dimensional magnets. With the short-time dynamic ap- 
proach, we determine the transition field H Cl and the static 
and dynamic critical exponents /?, v, z, £ and Q oc . The re- 
sults show that the fundamental Hamiltonian equation of 
motion belongs to a universality class very different from 
those effective equations of motion such as the QEW equa- 
tion and the Monte Carlo dynamics of the DRFIM model. 
In particular, the static exponent j3 of the </) 4 theory is 
much larger, and the global roughness exponent £ < 1. 
The local roughness exponent Q oc < 1 is compatible with 
experiments. 
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Fig. 1: The roughness function io 2 (t) is plotted at H=0.1 on a 
log- log scale. In the inset, the velocity v(t) is displayed on a 
semi-log scale. Dash lines represents constants. 
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Fig. 2: The velocity v(t, r) is displayed for different driving 
fields H on a log-log scale. For clarity, the curves of different 
lattice sizes at H c — 0.9828 are shifted down. The dash line 
shows a power-law behavior at H c . 
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Fig. 3: F(t) is plotted at H c = 0.9828 on a log-log scale. The 
dash line represents a power-law fit. 
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Fig. 4: The logarithmic derivative of v(t,r) is displayed at 
H c = 0.9828 on a log-log scale. The dash line shows a power- 
law behavior, and the solid line represents the fit with a power- 
law correction, d T In v ~ t 1 '" z (l + c/t). 
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Fig. 5: The roughness function Lj 2 (t) and pure roughness func- 
tion Dui 2 {t) are plotted against F(t) at H c — 0.9828 on a 
log-log scale. The dash line shows a power-law fit, while the 
solid line includes a power-law correction to scaling. 
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Fig. 7: The snapshot of the domain interface at H c = 0.9828. 
Black and white dots represent Sgn{4>i) — 1 and —1 respec- 
tively. 
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Fig. 6: The pure height correlation DC(r, i) is displayed for dif- 
ferent r at H c — 0.9828 on a log-log scale. Dash lines represent 
power-law fits. 
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